Criterion for explosive percolation transitions on complex networks 
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In a recent Letter, Friedman and Landsberg discussed the underlying mechanism of explosive 
phase transitions on complex networks [Phys. Rev. Lett. 103, 255701 (2009)]. This Brief Report 
presents a modest, though more insightful extension of their arguments. We discuss the implica- 
tions of their results on the cluster-size distribution and deduce that, under general conditions, the 
percolation transition will be explosive if the mean number of nodes per cluster diverges in the 
thermodynamic limit and prior to the transition threshold. In other words, if, upon increase of the 
network size n the amount of clusters in the network does not grow proportionally to n, the per- 
colation transition is explosive. Simulations and analytical calculations on various models support 
our findings. 

PACS numbers: 64.60.ah, 64.60.aq, 89.75.-k 
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The percolation problem, which deals with the struc- 
ture and connectivity of a network or lattice as some of 
its nodes or edges are removed, is one of the most widely 
studied problems in network theory P, During the 
last decade, much attention was devoted to the percola- 
tion transition on complex network, which has amongst 
others applications pertaining to virus spreading 0] and 
network failures [1, Q . 

Recently, Achlioptas et al. have discovered a new class 
of percolation transitions . The transition which they 
term explosive is marked by a vanishing "width" in the 
thermodynamic limit. According to the original Achliop- 
tas process, two candidate edges are selected at random 
at each timestcp. The edge resulting in the smallest prod- 
uct of its connecting cluster sizes is then effectively laid. 
Following the work of Achlioptas et al., explosive percola- 
tions have been observed on various networks and lattices 
using a wide variety of edge-addition rules, all of which 
involve the knowledge of the cluster sizes |7H13| . 

An important step towards a more profound under- 
standing of explosive transitions was made by Friedman 
and Landsberg Q. They put forward the criterion that 
a necessary condition for the appearance of an explosive 
transition is the existence of a nonzero fraction of nodes 
in the powder keg in the thermodynamic limit. In a net- 
work of size n, the powder keg F(t(n a ), n 1- ^) quantifies 
the number of nodes in clusters of size larger than n 1 ^ 13 
at time i(n CT ) when the first cluster of size n a makes its 
appearance. A trivial constraint is < 1 — j3 < a < 1 . In 
this Brief Report, we attempt to gain additional insight 
by presenting a more simplified criterion. Therefore, we 
consider the implications of a nonzero fraction of nodes 
in the powder keg on the cluster-size distribution for ex- 
plosive transitions. 

We assume that, close to the phase transition in the 
non-percolated regime, the number of clusters of size s, 
here denoted n s , is approximately described by a power 



law with a cut-off size s* for large s. Thus 

n, = As~ T e~ s / s * for s -» oo. 



(1) 



Here the normalization constant A is fixed by the con- 
straint n = sn s and, as occurs also for second-order 
percolation transitions [l[ , the cut-off length scales as a 
power of the network size; in other words s* oc with <; a 
non-negative exponent. Cluster-size distributions of the 
form of Eq. (JTJ) are for instance encountered for explosive 
transitions in simulations on the Achlioptas process on 
scale-free networks [§], on the model of D'Souza et al. [To[ 
and on the model of Cho et al. [HI [T3. as well as in the 
theoretical work of Da Costa et al. [14[- However, some 
studies reveal deviating cluster-size distributions |L5l - tl8| ; 
we will discuss them later on. 

Converting sums over the cluster-size distribution into 
integrals, the powder keg is given by 



F{t{n a ) 1 n 1 -f 3 ) 



sn, ds. 



(2) 



For large networks, the integral can be evaluated analyt- 
ically using the preceding assumptions and asymptotic 
expansions of the (incomplete) gamma function. We ob- 
tain: 
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F(t(n a ),n 1 ~ p ) oc n (1 ~^ )(2_r) if r > 2. (3) 



The fraction of nodes in the powder keg thus vanishes in 
the thermodynamic limit if r > 2. Conversely: 



F(t(n' T ),n 1 - f3 ) oc n if r < 2, 



(4) 
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in other words, the powder keg contains a nonzero frac- 
tion of nodes. In both regimes, the results arc verified 
numerically by an explicit evaluation of the integral in 
Eq. @. Note that the conclusions are independent of 
the value of <; and the exact definition of the upper and 
lower boundary of the powder keg, quantified by the ex- 
ponents a and (3. Therefore, we establish the criterion 
that the percolation transition is of the explosive type if 
t < 2 at the onset of the phase transition. 
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The actual determination of r in simulations is rather 
difficult. As an alternative, we present an equivalent cri- 
terion in terms of the mean number of nodes per cluster 
at the onset of the phase transition, (s) Q , the latter be- 
ing a more easily accessible quantity in simulations. Let 
n c be the number of clusters, then (s) = n/n c . Since 
n c = y^ „ n s , the mean cluster size (s) D can be evaluated 
analytically using Eq. ([T]). We deduce that in the ther- 
modynamic limit (s) diverges if and only if t < 2, a 
statement which can again be verified numerically. In 
other words, a percolation transition will be explosive if 
the mean number of nodes per cluster diverges at the 
onset of the phase transition, i.e. before the giant clus- 
ter is formed. For non-explosive transitions, a divergence 
of (s) can solely be caused by the formation of a giant 
cluster. 

Our criterion involving (s) quantifies the result of 
the intuitive mechanism underlying explosive processes 
known thus far: they all aim at avoiding any cluster 
size to become much larger than any other, therefore 
giving rise to a network with a small amount of clusters, 
all of which have more or less the same size [LlT - [l3j . 
Quantitatively, the transition will be explosive if the 
cluster density n c /n vanishes in the thermodynamic 
limit, i.e. if, upon increase of the network size n, the 
number of clusters does not grow proportional to n. 

We have checked our criterion, both numerically 
and analytically in various explosive models. Firstly, we 
examined the percolation process suggested by Cho et 



distribution follows approximately a power-law with a 
cut-off according to Eq. (p}. The starting point is a 
network consisting of n nodes without links. At each 
timestep, two clusters are selected, each with probability 
proportional to s Q+1 with s their cluster size and 
a G [—1,0]. The two selected clusters may be the same. 
A link is then laid between two random nodes of the 
selected clusters. Self-linking and multiple links between 
the same nodes are prohibited. The ensuing percolation 
transition was found to be exp losive if a < —0.5 and 
non-explosive otherwise [HI, [Hj. Fig. [T] shows the mean 
number of nodes in a cluster at the appearance of the 
first cluster of size ^/n, which is generally considered as 
the onset of the phase transition. 

The non-explosive regime with a — corresponds 
to random percolation. An expression for the cluster- 
size distribution for random percolation was obtained in 
Ref. fl9j . An exact computation of the mean cluster size 
per cluster in the thermodynamic limit then reveals the 
finite asymptotic limit (s) Q = 2, which nicely corresponds 
with our simulation data as seen by the dots and dotted 
line in Fig. [1] 

Also in the explosive regime, a = — 1, the cluster-size 
distribution can be obtained analytically [ill ]. At time t, 
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FIG. 1. (Color online) Log-log plot of (s) , the mean num- 
ber of nodes per cluster at the onset of the phase transition, 
as a function of the network size for various percolation pro- 
cesses. The onset of the phase transition coincides with the 
time when the first cluster of size appears. For the Cho- 
Manna process (CM), results for an explosive regime [a, = — 1, 
green right triangles) and a non-explosive regime (a = 0, 
black circles) are shown. The (black) dotted line shows the 
analytically obtained asymptotic limit for the non-explosive 
regime, while the (magenta) dot-dashed line gives the analyti- 
cal result for the explosive regime. Also shown are simulation 
results for the explosive min-cluster sum rule (MC-10, blue 
up triangles) and the explosive extended-product rule (PR- 
10, red down triangles). All data points are averages over 20 
realizations of the percolation process. 



al. 11( and Manna et al. 12j, in which the cluster-size anc j thus 



1 - t 



(6) 



An expression for the mean number of nodes per cluster 
at the onset of the phase transition (i.e. when there is 
exactly one cluster size larger than n 1_/3 ), is then found 
by solving the equation 



1= £ n *(*) 



(7) 



n s (t) = n(l - t) 2 t 



2 + s-l 



(5) 



for t and introducing the result in Eq. (|6|). Results of 
a numerical evaluation of this procedure for j3 = 0.5 
are shown in Fig. [TJ which also shows the value of 
(s) determined in the simulations. The theory is 
found to agree well with the simulations. Both reveal 
a mean cluster size per cluster which diverges upon 
approach of the thermodynamic limit as (s) Q cx n 1 
with 7 = 0.42 ± 0.02. We conclude that simulations 
and analytical computations on the Cho-Manna model 
confirm our second criterion. 

Our presented analysis and simulations so far re- 
lied on the ansatz for the cluster-size distribution, 
Eq. (Q]). Based on the cluster-size distribution near 
criticality, different types of models exhibiting explo- 
sive percolation may be distinguished. A power-law 
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distribution with a bump is observed in the models of 
Refs. [HI, [l6|, while a lattice model yields a Gaussian 
cluster-size distribution [171 ] . On the other hand, there 
exists an extremely explosive model in which the cluster 
size distribution is not continuous but highly discrete. 
In the latter model, at each timestep maximally two 
cluster-sizes are present in the network fl8j . 

We have performed simulations on two models in which 
the cluster-size distribution features a bump for large 
cluster sizes. Both models start from n distinct nodes. 
At each timestip, ten possible edges are selected at ran- 
dom. In the first model, a straightforward extension of 
the rule of Achlioptas denoted as PR-10, the edge which 
results in the smallest product of its connecting clus- 
ter sizes is effectively laid. In the min-cluster sum rule 
(MC-10), out of ten edges, the edge which minimizes the 
size of the component formed if the edge is occupied, is 
laid. The cluster-size distribution in the non-percolating 
regime does not follow exactly Eq. ([1]) [la] : although 
it is characterized by a power-law behavior for small 
cluster sizes and an exponential cut-off for large sizes, 
there is a pronounced bump in the probability distribu- 
tion for intermediate sizes. Nevertheless, as shown in 
Fig. [TJ for both processes the mean number of nodes per 
cluster follows a diverging power-law of the network size, 
(s) cx rf. For the MC-10 rule, we obtain 7 = 0.44±0.02; 
while 7 = 0.47 ± 0.01 for the extended product rule. Al- 
though the derivation presented in this Report is based 
on Eq. (fT]), it seems that our criterion holds in explo- 
sive models in which the ansatz is no longer true. In our 
opinion, the models with a bumped cluster-size distri- 
bution still satisfy our criterion since their distributions 
contain two essential ingredients of distribution Eq. ([1}, 
that is, a power-law distribution at low cluster-size and 
an exponential cut-off length which diverges in the ther- 
modynamic limit. 

We have also tested our criterion using a model with a 
highly discrete cluster-size distribution, based on an ex- 
treme form of the Achlioptas process. The PR-00 rule 



takes into account all possible links and then selects the 
link which results in the smallest product of its connect- 
ing cluster sizes. Suppose for simplicity the system size n 
to be an integer power of two. The cluster-size distribu- 
tion is then trivial since maximally two types of clusters 
exist at each timestep [lH . Suppose that 2^ is the largest 
power of two smaller than n 1- ' 3 , thus 2^ 
At the moment of the first formation of a cluster of size 
n i-p or larger, there is a single cluster of size 2^ +1 , while 
all other nodes belong to clusters of size 2*. Since there 
are 2~^n — 2 such clusters, we obtain 



In the thermodynamic limit, the mean number of nodes 
per cluster diverges at the onset of the phase transition. 
This exact calculation thus shows that our criterion is 
also valid in a model in which only a finite number of 
cluster sizes are present. 

In sum, under general conditions a percolation process 
will be explosive if the mean number of nodes per 
cluster diverges at the onset of the phase transition in 
the thermodynamic limit. Equivalcntly, if the cluster 
density n c /n vanishes in the thermodynamic limit prior 
to the critical point, the transition is explosive. Al- 
though the criterion is deduced with a specific ansatz for 
the cluster-size distribution, simulations and analytical 
calculations on various explosive models indicate a more 
general validity. 
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